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^ ! Abstract 
& ■ 

The complete form of the high-temperature expansion of the one-loop contribution to the free energy of a 
scalar field on a stationary gravitational background is derived. The explicit expressions for the divergent and 
finite parts of the high-temperature expansion in a three-dimensional space without boundaries are obtained. 
These formulas generalize the known one for the stationary spacetimc. In particular, we confirm that for 
a massless conformal scalar field the leading correction to the Planck law proportional to the temperature 

Ih ' squared turns out to be nonzero due to non-static nature of the metric. The explicit expression for the 

bJQ' ..... , 

so-called energy-time anomaly is found. The interrelation between this anomaly and the conformal (trace) 

t^J" ' anomaly is established. The natural simplest Lagrangian for the "Killing vector field" is given. 
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I. INTRODUCTION 



The high-temperature expansions of the partition functions are the classical subject of quantum 
field theory on a curved background. As for the one-loop contributions to the free energy of 
quantum fields, we just mention the works [[j|-|lO|| where the different approaches to this problem were 
implemented. Surprisingly, in spite of the fact that the first attempts to find the high-temperature 
expansion were undertaken over about thirty years ago Q, the probl em was not completely solved, 
to our knowledge, in its general statement for the arbitrary stationary (non-static) gravitational 
background. The present paper is aimed to fill this gap. Such a high-temperature expansion for the 
free energy will be derived here. In particular, for a massless conformal scalar field we shall obtain 
the leading correction to the Planck law proportional to the temperature squared and confirm the 
result of H |i~0| . This correction is absent for a static gravitational background [ffl, ^, Also we 
shall derive the explicit expression for the finite part of the high-temperature expansion, which 
seems to be a new result. 

Apart from the immediate implications for the Casimir effect, astrophysics, and cosmology, the 
high-temperature expansion of the free energy for fermions can be used to analyze the derivative 



expansion of the one-loop contribution to the effective action at zero temperature (see, e.g., [11, ft 



regularized by the energy cutoff. According to the general prescriptions of the renormalization 



theory [13|, the structures appearing as the divergencies in the effective action must be included to 
the initial action to cancel these divergencies (not necessarily to zero). As a rule, these divergencies 
and the finite part depend nontrivially on the Killing vector defining the stationarity of the 
background and the vacuum state of the quantum fields (||14]|, see also [15-17]]). We shall see in the 
present paper that this is indeed the case. Hence, the analysis of the high-temperature expansion 
sheds a light on the low energy quantum dynamics of the vector field £^ as it was discussed in 
p2|, |14|. To see how this works, one may bear in mind the expansion of the Heisenberg-Euler 



effective action [18, 19 1, where the first non-trivial term describes the light by light scattering. 

Another interesting point following from the results of the present paper is the interrelation 
between the conformal (trace) anomaly and the energy-time anomaly (the notion of the latter was 
introduced in ||14j| , see also below). As we shall see, we cannot renormalize the quantum theory of 
a massless conformal scalar field in such a way that both the conformal and energy-time anomalies 
vanish. The elimination of the one anomaly results in the appearance of the other and vice versa. 
It is noteworthy that the factor at the logarithm of the temperature entering the high-temperature 
expansion and determining the anomaly turns out to be independent of the Killing vector and 



coincides with the standard expression for the conformal anomaly ||, 1C|. However, the finite part 
of the high-temperature expansion does involve certain contractions of the Killing vector field. 
Loosely speaking, the problem of dependence of the effective action on the Killing vector comes 
from the infrared modes of the gravitational field and so it is reasonable that the renormalization 
group beta function defining the ultraviolet behaviour of the theory is independent of the Killing 
vector. Though, of course, this fact is not a priori obvious since both in the ultraviolet and the 
infrared regions we have the mode functions defined with respect to the same Killing vector (see, 
however, [ |lO|| for the possible proof). 

The paper is organized as follows. Since the subject of the article is rather technical, we par- 
tition it into the small sections distinguishing the major successive steps of the derivation, which 



culminates in Sec. VII where the results and its implications are discussed. In Sec. |IT|, we provide 



general formulas for the one-loop correction to the free energy. The main technical tool, which we 
shall employ to derive the high-temperature expansion, is the heat kernel of a Laplacian operator 
(see for review). Therefore, in Sees. Ill, IV, and ^ we reduce our problem to the evaluation of 



the heat kernel expansion coefficients. At first, Sec. |T|, we reduce the three-dimensional problem to 
the four-dimensional one to provide the explicit general covariance to the expansion. The method, 
which is used here, is a finite dimensional analog of the gauge fixing procedure in the functional 



integral. Then, in Sec. [V, we resum the heat kernel expansion applying the theorem proved in [21]. 
Notice that we do not use the conformal transformation to derive the high-temperature expansion 
H, fl|-[ro|l, but apply the method directly to the Fourier transformed Klein-Gordon operator. So, 
our approach is rather close to the one used in g]. After that, the problem becomes in essence 
the same as in a flat spacetime. In Sec. |V|, we evaluate the expansions of necessary integrals em- 



ploying the procedure used in [11]. Then, in Sec. VI, a general formula for the high-temperature 



expansion of the free energy is obtained. In Sec. VII, we particularize the general formulas to 



the three-dimensional space and single out the divergent and finite parts of the high-temperature 



expansion. At this point we essentially employ the results of [22|, where the heat kernel coefficients 
were derived up to a§ . In Appendix |A|, we give some useful relations for the metric possessing the 
Killing vector. In Appendix [B], the relevant parts of the heat kernel coefficients borrowed from 
[20|, H-|H are calculated. 

Despite the fact that we study the high-temperature expansion of the free energy of quantum 
scalar field on a gravitational background, the results are easily generalized to higher spins and 
other stationary backgrounds. The problem lies only in the amount of calculation, which increases 
due to the additional terms in the heat kernel expansion. A non-zero chemical potential can be 



also included. We postpone the investigation of these problems for the future research. 

We shall use the following conventions for the curvatures and other structures appearing in the 
heat kernel expansion: 

R a piiv = ®\pFu]P + ^7^3]^' = Oilcan R = Rp- W 

The square and round brackets at a pair of indices denote antisymmetrization and symmetrization 
without 1/2, respectively. The Greek indices are raised and lowered by the metric g^ u which has 
the signature —2. Also we assume that the metric possesses the timelike Killing vector 

jo^u = o, e = g^ee > o. (2) 

The space dimension will be denoted by d and d = 3. Nevertheless, we shall put d = 3 only in the 
final result. The system of units is chosen such that c = h = 1. 



II. GENERAL FORMULAS 



Consider a scalar quantum field on a stationary gravitational background at a finite reciprocal 
temperature /3. The free energy for this system is defined in the standard way 

e^ F :=Tve^ n , (3) 

where % is the Hamiltonian of the scalar field expressed in terms of the creation-annihilation 
operators. The mode functions of this field corresponding to the energy u span the kernel of the 
Klein-Gordon operator, 

\g\ l ^{x)\g\ L /*{y) 

= \g\~ l/ \x) 

where all the time derivatives should be replaced by — iu. This operator, which we denote as H(ou), 
must be supplemented by the appropriate boundary conditions. To simplify further calculations 
we assume that the system considered is large enough to neglect the boundary effects or the space 
represents a compact manifold without boundary. The operator H{oj) is Hermitian with respect to 
the measure y/\g\ on the square-integrable functions depending on x. 

Our aim is to calculate the one-loop contribution to the free energy (||). To this end, we use the 
fact that 

d u Tt d 9(H(u))) =J2^{ef k (u))6(u-u h ) (5) 
k 
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-1/4 



(x)du 



\9\g 



-1/4, 



X) 
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K x - y) 
\g\ 1/4 (y)' 



(4) 



defines the signed spectral density in the w-space for the self-adjoint operator H(uS) possessing the 
eigenvalues £fc(w). Here the Heaviside step function 8(H(u) — e) is the spectral decomposition of 
unity associated with H(co) and ui^ is the solution to the equation 

e*(w) = 0. (6) 

If this equation has several solutions then the sum over all such solutions should be taken in (||). 
Note that u may be any parameter of the self-adjoint operator and not only the energy. In the case 
when uj is the energy as defined above, s' k (uj) is usually positive for to > (the particle branch) and 
negative for oj < (the antiparticle branch). Differentiating the equations 

H(u> k ,m 2 )ip k (m 2 ) = 0, (V>A#fc) = 1, (7) 

with respect to m 2 , it is easy to see that the last statement is valid for the same cases, when 
the standard prescription m 2 — > m 2 — iO defines the Feynman propagator. It follows from the 
quasiclassical asymptotic of the operator H(u), where all the derivatives entering it are replaced by 
the momenta —ip^, that E k (oj) can be negative for to > under the ergosphere, i.e., the antiparticle 
branch appears in the positive-frequency region. This results in the particle creation p5[ and the 
reconstruction of the vacuum state. Therefore we restrict our consideration by the region out of 
the ergosphere. 

Assuming that the lowest particle energy is strictly positive (this is the case when the system is 
placed in a finite "box"), we can write the one- loop contribution to the free energy from the particles 
as 

T f3F= / 00 ^Tr^(FH)ln(l±e-^) = ±/3 f°° du P*^}"^ , (8) 
jo Jo e ± 1 

where plus corresponds to fermions and minus is for bosons. The contribution from the antiparticles 

is given by the same integral, but with H(—lo) instead of H{uS). Hence, the one-loop contribution 



to the effective action from one bosonic mode at zero temperature reads (see, e.g., [11, |12|| ) 



/3-K) 







/3-s>0 



where plays the role of the energy cutoff and T is the time interval tending to infinity. The 
fermionic free energy for a scalar field (i.e. (g) with the plus sign) can be also used to estimate the 
partition function for fermions when the spin-gravity interaction is negligible. For further purposes 
it is useful to represent the step function in the integral form: 

9 <"H) = / o 5-— . do) 



where the contour C runs along the imaginary axis from top to bottom and encircles the origin 
from the left. The exponent on the right-hand side is the so-called heat kernel. In order to obtain 
the high-temperature expansion (/? — > 0), it will be sufficient to know the expansion of the trace of 
the heat kernel in a series in r. 



III. REDUCTION TO d+l 

Of course, we can apply the formulas from 22-p4|| directly to (|To|), but this leads to rather 
cumbersome expressions that are not explicitly covariant because of the (3 + 1) decomposition. 
Therefore, we reduce our 3-dimensional problem to the 4-dimensional one preserving explicitly the 
general covariance. To this end, we introduce the auxiliary integration variable, 

e-^W = f°° dp e- rH ^8(p - u) = f°° ^ e" T ^ (po)+ ^ (aj2 " P °%(po - (11) 

J —oo J —oo 

and employ the standard Gaussian representation of the delta function 

S( Po -u))= lira W- e -^(po-o,) 2 r i/ 2 . / 12 j 
A— ¥-oo V vr 

The second term in the exponent in ( pi] ) tends to zero when A goes to minus infinity. The usefulness 
of such an addition will become clear soon. The expressions with a finite A correspond to the same 
system but with a Gaussian broadening of the spectrum. Denoting 

H x ( Po ,oj) :=H(p ) + X( Po -uj) 2 + ^(uj 2 -p 2 ), (13) 

we have 

Tr d e~ TH{uj) = lim J- f°° d Po Tr d e -^A(Po^) T l/2 = 1 Um 2v ^A Tr D e^^r 1 / 2 , (14) 

A— >-oo V VT J^oo T A— ¥— oo 

where D := d + 1 is the spacetime dimension. Then we make a similarity transform such that the 
diagonal of the heat operator in the x-representation and, consequently, its the trace are left intact, 
while the Hamiltonian becomes (see also Appendix |A| for the notation) 

H G := {ef ,A H x (er 1,A = -GT{V» - i^){V v - iug„) - X, (15) 

where := ^/S, 2 and 

G„ v = 9fip - e^9v + A" GT = gT - £ 2 (1 - A£ 2 ) 5 V, det = X^g 2 det g^, 

X = V-E, F = i^_iv%, £=^-m 2 , 

(16) 

and hn := In y^ 2 "- The connection V M is the Levi-Civita connection compatible with the metric 
G^y. The latter is negative definite so that the operator Hq is of the Laplacian type. 



IV. RESUMMATION OF THE HEAT KERNEL EXPANSION 



Now we can use the heat kernel expansion 

oo 

(4vr) D / 2 ' 



oo k—D/2 



k=0 

It is also relevant here that the space has no boundary. Otherwise the terms with the half-integer 
k appear. To exploit the heat kernel expansion, we resum it using the theorem proved in pi] ] 

oo k—D/2 

( X \e- TH 0\y) = e r[X( X )^R(x)] {xle -r[X( X ^R( X )] e -rH G \ y) = & t\X (x)+$R(x)] y ^ fa y) ) 

fo ( 4vr ) ' 

(18) 

where dk are obtained from a k by removing the terms containing X and R without derivatives. As 
for the X terms, this theorem simply follows from the observation that 

fa e -T[X(x)+±R(x)] e -TH G ^ ( 19 ) 

is invariant under the transform X — > X + const and so its expansion does not contain X without 
derivatives. On the other hand, expanding 

(x\e- T l x M + *^e- TH °\y) = e -^W+^Wl {x\e^\y) (20) 

in t, we see that at the same power of r: 

a k = a k + 0(X,R), (21) 

where 0(X, R) are the terms containing X and R without derivatives and a s with s < k. If dk does 
not contain X and R without derivatives in virtue of summation of such terms to the exponent then 
dk is obtained by a mere obliteration of the terms in a k proportional to X or R without derivatives. 
The fact that dk do not depend on the scalar curvature R can be proven by a direct inspection of 
the heat kernel expansion coefficients [^] or by solving the defining equations for the heat kernel 



in the Gaussian approximation (see, e.g., ]14|). 



V. EVALUATION OF THE INTEGRALS OVER r AND u 



For brevity, we introduce the notation 



rh 2 := m 2 + V + \k = m 2 + -h 2 + \{R - V%) - ^(A" 1 - f)f 2 = 

24 ^2) 

= m 2 + ^2 + 1 [R _ (A -1 _ f)R^g- - A-VV%], 



7 



where we have used the relations from Appendix [A]. The quasiclassical expansion for the free energy 
@ after the above transformations becomes 

- F = ^ lim JEf* [ dx^G r [ j^T^-^V-^- 1 , (23) 

where dk{x) := dk(x,x). Of course, on using the heat kernel in the form (|l~8|), we neglect all the 
exponentially suppressed at /3 — )■ terms in the high-temperature expansion. For the complete 



analytical structure of the heat kernel in the r-plane see, e.g., [14|. The integral over r can be 
easily taken 

a We 1 = vm-k + iA e ~ m ) ^ " m ) ■ (24) 

Comparing (p4[), (p3[), and (Jl6|) , we see that the imaginary units are all canceled out and the 
expression ( p3| ) is real as it should be. The formula (p4|) ought to be understood in a distributional 
sense for k > d/2. The simplest way to take this into account is to assume that the variable d is 
complex and tends to its physical value 3. This rule can be check by convolving ( p4[ ) with a test 
function single-valued near the real positive semiaxis u > 0. Only for those test functions does the 
prescription of the analytical continuation work. In our case this requirement is evidently fulfilled. 

Now the integral over u has to be evaluated. The heat kernel expansion coefficients 5fe(o;) are 
polynomial in oj: 

2k 

a k (cj,x)=:J2^ ] ^)(9 2 ^) j/2 , (25) 

3=0 

where the upper summation limit is dictated by dimensional reasons. Observe that the heat kernel 

(x\e- THG ^\y) (26) 

is Hermitian at the real r and therefore its diagonal is real in this case. On the other hand, from 
Eq. ( |l~5l) we see that the complex conjugation effectively results in a change of the sign of ui. Hence, 
by the uniqueness of the analytic continuation in r, the diagonal of (p6|) is an even function of ui 
for any r. Then the polynomial ( p5|) contains solely the even powers of the energy oj as long as E 
is an even function of oj. 

So, we need to take the integral 



I 



Bdw 6{oj 2 g 2 - fh?) , 9 n xr//9 i..,-/?/, rh 2 \d/2-k 

et 1 1 TW- k + 1) (g " ' 7 "V ~ Z&) ■ P7) 

The further procedure is quite analogous to the case of a flat spacetime (see, e.g., 26-28]]). We 
shall follow Substituting the expansion 

m 2 \ d/2—k _ ~ (_i)"r(d/2 - k + 1) / m 2 y 

0J 2 g 2) 2^ Q nT(d/2-k-n + l)\oj 2 g 2 ) [ j 



to the integral /, we reduce it to a sum of the incomplete zeta functions. Furthermore, we expand 
these incomplete zeta functions making use of the formulas fl29|| : 



doj 



&" — 1 



{-l) n C{-n)a u+n 
x T(n + l)(u + n) ' 



The term with n 



« (-l)"m d - 2fc V 



n=0 



^— = (1 - 2 1 -)I»<» - V (1 - 2 1+ »)ffi^. 
e w + l v ' usu ^ v J T(n + l )(v + n) 

n=0 

1 is understood as the limit n — > — 1. For the bosonic case we obtain 

r(d + j - 2jfc - 2n + l)C(d + j - 2k - 2n + l){f3 T m) 2k+2n ~ d 

(-l)'CH)0W +t+1 



(29) 



n!r(d/2- fc-n + 1) 



E 



. (30) 



r(Z + l)(d + i - 2k - 2n + I + 1) 

where /3t '■= \f£?P is the Tolman reciprocal temperature and we have assumed that rh 2 > 0. The 
sum over n of the second term in the square brackets can be taken and is expressed through the 
beta function: 



•i 



n=0 

Whence we get 



n\T(d/2 - k - n + \){d + j - 2k - 2n + 1 + 1) 



k- {d + j + l + l)/2) 

2r((i-j-z)/2) 



(31) 



-l)"r(d + j - 2fc - 2n + l)C(d + j - 2k - 2n + 1) . 2n „ 2 k+2n- j -d 



+ 



n=0 



n!r(d/2 — k — n + 1) 

, (-i)'C(-Qr(fc - (d + j + j + 1)/2) w?+m flt+1 
^ zr(i + i)r((i-j-i)/2) 

So, we have all what we need to obtain the high-temperature expansion. 



(32) 



VI. HIGH-TEMPERATURE EXPANSION 



Substituting the expression (32) to (|23|), we eventually arrive at 

- °° (-l) n+k T(d + j-2k-2n + l)Q(d + j - 2k - 2n + 1) 



k,j=0 J n=0 



W(d/2 -k-n+ l)(3p j - 2k - 2n+1 m- 2n 



+ 



+ £ 



n 



(-i)'+ fc C(-Qr(fc - (d + j + 1 + 1)/2) d+j _ 2fc+w 



K=-l 



2r(z + i)r((i-i-/)/2) 



/St » (33) 



where the limit of the infinite A is implied. The expression standing in the square brackets is regular, 
when d tends to a positive integer, since ( |30| ) is the entire function of d. Therefore, we can set d 



9 



to its physical value in the factor at this square brackets. Besides, due to the property of the zeta 
function a "half" of the terms in the sum over I vanish. These are the terms with I = 2, 4, 6, . . . As 
far as the fermions are concerned, similar calculations lead to (to 2 > 0) 

oo . 

- F f = E j dx 



( 4vr )rf/2 



fc,i=Q • 

x r Vn 92n + 2fc-j-^ (-l) n+k T(d + j-2k-2n + l)C(d + j - 2fc - 2n + 1) 
^ j n!r(d/2-fe-n + l)/3^- 2fe - 2 ^ 



fn ni+^ (-i) i+fc C(-/)r(fc - (d + j + 1 + 1)/2) d+? -_ 2fc+ , +1 , 
^ ( } 2ra + i)r((i-,-/)/2) m ^ 



(34) 



The terms at the fixed power of f3 both for bosons and fermions are the same up to an overall 
numeric factor and in that sense are universal. However, the expansion of the free energy for 
bosons involves one additional contribution with I = — 1. 

Let us analyze the first and second terms in the square brackets of ( |33"|) and (34) separately. We 



see that both for bosons and fermions the second terms in the square brackets are not expanded in 
the increasing powers of /?. Rather, we have a derivative (or large mass) expansion of the terms at 
the every fixed power of the reciprocal temperature. It is not difficult to write out a closed form 
expression for a whole series at the fixed power of f3. With this end in view, we introduce the 
function (cf. § 0) 



al(m 2 ) 



[ dssU r dLoJTr d e- sH ^. (35) 



The function under the integral over s given by the integral over u is analytic for Re s > 0. It should 
be analytically continued to the imaginary axis and to the vicinity of the origin where the contour 
C lies. Usually, this can be achieved by rotating the integration contour in the w-plane (see the 
high frequency asymptotics in Eq. (^3|)). Notice that the function ( |35| ) is not the generalized zeta 
function as long as H(uj) possesses the negative eigenvalues. Although these functions are closely 
related. Repeating all the above calculations, one can convince oneself that the expansion (^) can 
be rewritten as 

-F- V ! Ac (-l)*r(d + j - 2k - 2n + l)C(d + j - 2k - 2n + 1) 

" k f^ = J X (4vr)<V2 n W(d/2-k-n + l)4 +j - 2k - 2n+1 (-m*)-" 

+ E r(z + i) (36) 

where e = (d — d)/2 is the complex number tending to zero and d is the physical dimension {d = 3). 
The analogous formula holds for fermions as well. The function a l v {m?) depends only on to 2 and 

10 



the background fields. Consequently, we need certain additional assumptions apart from j3 — )■ 
in order to obtain an explicit expression for it. Henceforth we assume that the gravitational field 
varies slowly such that the derivative expansions presented in (^) and (Q) make sense and provide 
reliable approximations for the functions (|35|). This holds when the characteristic length of the 
variation of the gravitational field is much larger than the Compton wavelength associated with 
the effective mass m. In any case, these second terms give subleading contributions to the high- 
temperature expansion. For fermions, they stand at the nonegative integer powers of f3, while for 
bosons the leading contribution from these terms is proportional to /3 . The first terms in the 
square brackets in the expansions ( |33[ ) and (|34|) are much more singular at (3 — > 0. 

Consider the first terms in the square brackets of Eqs. (33) and (p4|). There is a finite number 



of such terms at any fixed power of (3. Indeed, at any fixed number k, the number j must be less 
or equal to 2k by dimensional reasons. However, inasmuch as we resummed the expansion and 
collected all the terms without derivatives of E to the exponent, the latter terms are absent in a^. 
Consequently, the worst terms, which contain the maximal power of lo at the fixed dimension 2k, 
are of the form 

VEVE ■ ■ ■ VE. (37) 



If the dimension 2k is not a multiple of 3 then the worst terms look like ( p7[ ) , but with the additional 
covariant derivative, or with f2 2 (see, for example, ( B17|) ). Thus, we conclude that (see also g, |30| |) 



j < [4fc/3], (38) 
where the square brackets denote the integer part of the number. These worst terms stand at 

and so the power of (3 increases. Further, we restrict ourself by the singular and finite parts of the 
high-temperature expansion as they give the leading contribution to the free energy. From (|3^) we 
conclude that the maximal number k that we need is determined by the equality 

2k - [4k/3] = D, (40) 

and for a three dimensional space k < 6. Fortunately, the heat kernel coefficients with k < 6 
were found in |22|| and we just can borrow these results. The relevant parts of these coefficients are 
presented in Appendix [B[ 



Now we turn again to the second terms in the square brackets of the expansions (J33|) and (34). 



The contributions with the positive even / vanish. Moreover, the contribution from these terms is 
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zero when I is an odd positive number and j is an even nonegative number. As we discussed above, 
the number j must be even and nonegative. The explicit expressions for a& given in Appendix |B| 
are found to be even functions of uj for those which we are interested in and confirm thereby the 
general considerations. 

Note that these second terms were disregarded in [^], ^, [TIJ . It is relevant for a correct evaluation 
of these terms that we did not make a conformal transformation to the optic metric and resummed 
the heat kernel expansion. The resummation gives a non-perturbative expression for the heat kernel 
which is more reliable in terms of the original metric rather than the optic one (see for details p4]|). 
In principle, the first contributions in the square brackets of Eqs. (33) and ( |34] ) may be evaluated 
in terms of the conformally transformed metric. However, in order to provide a transparent way 
to the poles cancelation in the d-plane they have to be expressed through the same metric as the 
second terms. 



VII. RESULTS AND IMPLICATIONS 



Bearing in mind these observations, we can write for d = 3 in the bosonic case 



dx.yj\g\ ^ k r w T(D + j - 2k - 2n)((D + j - 2k - 2n)d 



2\-n 



^ d/2 k% 1 n\T{d/2 -k-n+ l)$ +j - 2k - 2n (-m 2 ) 
(_ fh 2\{D+j)/2-k~i3) m 2 Ble 2 ~< 

+ w^hwmhm [ ln ^ - ^ + >v* - * + d + - m 

+ 5, T(k - d/2)^m d - 2k - a?m D+ ^ 2k T{k ~} D \ j ) /2) }, (41) 
J,U 1 1 J 2P T k 4r((l-j)/2) V v ' 

where 7 is the Euler constant. The prime at the sum over n says that the singular terms are 
discarded, the second term is zero by definition when the argument of the factorial is negative, 
and the last term vanishes by definition when the gamma function entering the numerator tends 
to infinity. The last term stands at the negative power of rh. Later on we shall cast out such 
contributions although it is these contributions which are "protected" from the high energy physics, 
i.e., the particles with a large mass (but with m -C /3 X ) give a negligible contribution to these 
terms. The high-temperature expansion in the fermionic case looks like 



-F* 



[ dxVW\ A k r 2 2k+2n-i- d) T (D +3 ~2k- 2n)gP + j - 2k - 2n)d 

J M d/2 k % [ ' ' n\T(d/2 -k-n + l)^- 2k - 2n (-rh 2 )- 



^2)(D+i)/2-fc a Cf) Am 2 p 2 e 2 ^ 



-m 



4((D+j)/2-k)T((l-j)/2)t tt 
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[ In H(D + j)/2 - k + 1) + m - j)/2)] 



+ k 4r((i-i)/2) /• (42) 

The digamma functions appearing in the expressions are easily calculated using the formulas [ pEj] : 

1 n 2 n 1 

^2~ n ) =E^3T- ln4 - 7 ' ^(l + n) = ^-- 7 . (43) 
fc=l fe=i 

More specifically, taking into account that a® = 1 and a\ = 0, we derive for bosons 

-f„ = / ^-M{^f - §jW - ¥<? - W) + f&(™ 3 + + ■ • •) 



3 



1287T 2 

and for fermions 



[m 4 - m\l~af + §#) - §(# + 2af + f af + 223®)] +...}, 



In [*< + ^ + 14 4) ) + 2af + + 14 4) + ¥4 6) + ¥ 4 8) ] (45) 

3 f^4 ™2^2~(2) . 5 ~(4)n 4/~(2) , ~(4) . 23~(6) . OQ -(8), 



. [rf _ ^2 ( | 5 (2) + | 5 (4) } _ 4 (5 (2) + 2 ~(4) + ^(6) + ^-(S))] + j 



1287T 2 

The dots in the expressions above denote the terms at the negative powers of rh or the positive 
powers of /3. The coefficients diverging at (3 — > provide the energy-time anomaly ([|l4]] and also 
[ p~6| - |l7| |) when the energy cutoff is used for the regularization of the one-loop contribution to the 
effective action. Recall that the energy-time anomaly represents the variance of the renormalized 
effective action under the uniform dilatations of the Killing vector (see below) . The flat spacetime 



limit of the expansions ( |44|) and ( |45|) coincides with the known one |26-28|. 



Now we consider the coefficients at the different powers of (3 in detail. First, we observe that 
j <&vWV - \af - \af) = /T 2 / d X ^\[(m 2 + \ R )C 2 + - ^(g%)]. (46) 

If we omit the last term representing a total divergence then we confirm the results of |l], ^, || that 
the term at /3 -2 is absent for a massless conformal scalar field on a static gravitational background. 
Meanwhile, we see that this is not the case for a stationary spacetime. For such spacetimes, the 
leading correction to the Planck law is of the order of f3~ 2 and proportional to the Maxwell-like 
term f 2 . Formula (46) coincides with Eq. (4.18) of Q up to the integral of a total derivative. 



Under certain circumstances the last term in (46) reducing to a boundary term also considerably 



contributes. However, in this case other boundary contributions that we discarded earlier must be 
included into the free energy too. 
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In the bosonic case, the next correction is proportional to For a massless scalar field on a 

curved background satisfying the vacuum Einstein equations K^ v = 0, this contribution becomes 
purely imaginary as long as the correction to m 2 is negative in this case |14|. There is also an 
imaginary contribution to the free energy due to the logarithmic term. The imaginary contributions 
to the effective action signalize about a certain instability of the vacuum of quantum fields which 
was initially chosen. As a result, we may anticipate that a non-zero vacuum expectation value of 
the scalar field develops (the Bose-Einstein condensation). 

Of course, if we accurately evaluate the high-temperature expansion of (||) then there are not 
any imaginary contributions to the free energy as the expression (H) is real-valued. For fh 2 < 0, the 
additional terms appear in (|D and (fH such that the imaginary part of the free energy vanishes. 



In the fermionic case, the contributions to the real part of the free energy standing at /3 -4 , (3~ 2 , /3°, 
and ln/3 do not change, i.e., the additional terms cancel exactly the imaginary part coming from 
the logarithm and add certain contributions at the higher powers of f3. At the same time, in the 
bosonic case, the imaginary part is also canceled out, but the infrared divergence arises in the term 
with /3 . Other contributions to the singular and finite parts of the high-temperature expansion 
remain unchanged. The term at entering the free energy Fj, becomes 

00 l 1\k ( \^\d-2k%(°) \2 



r /u ( -\\ k ( \™\ d - 2k n vu > \ 2 

I E m-2 - a i) i^F + 7 + m2 - k+1)] 

\m\ d+j - 2k T{j/2) ] f (ixy^ 
T((d + j)/2-k + l) J J 12n 2 fr 



k=0 

[ ST ] ~U) ( lV /2 \rh\ d+j ~ 2k T(j/2) } f^yfr |3 . g 2 X 2 e^ 3\rh\ (2) _ (4) , 

g 4 r(fd+.7)/2-fe+i) f = ./ i2^H ln ^r + — {a2 + °3 )+ - 



(47) 

where A is the infrared energy cutoff. The contributions vanishing at A — > were neglected and 
dots denote the negative powers of fh. 

These corrections to (44) and (45) have their origin in the fact that at fh 2 < the integral (p] 



can be represented as the sum of two integrals along the contours [A, iy^|m|] and [iv^^l^l' +°°)- 
The latter integral was already evaluated and its contribution to the high-temperature expansion 
is given by ([|l]) or (^2|). The former integral, 

rWPW p duj _ 7h 2 ) d/2-k 
h e#"±l r(d/2-ifc + l) ' ( } 

is readily taken if one expands (e^ u =t 1) _1 in a Laurent series in fiuj. The infrared cutoff is only 
needed in the bosonic case for the contribution at (/3w) _1 and j = 0: 

p duj ^ g 2J2 _ ^d/2-k riy/e\m\ ^ ^2J2 _ ^d/2-k 



J- x ePu-1 T(d/2-k + l) ~ J~ x u T(d/2-k + l) 
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(49) 



This integral, in turn, can be estimated at A — > making use of the expansion of the incomplete 
beta function p9| : 

n=0 (50) 
J a x rail [p — n)n 

To sum up, we see that at m 2 < the infrared divergencies arise in the free energy. Therefore, 
a more accurate analysis of the low energy modes is required. These divergencies may be cured by 



introducing the Bose-Einstein condensate and summing the ring diagrams (see, e.g., [26, 28]) or, 
as in the flat spacetime, by introducing the boundaries and accounting for the finite size effects. In 
the case, when 

/ < m- 1 « L, (51) 

where I is the characteristic scale of the variations of the gravitational field and L is the size of 
the "box", the boundary effects are negligible and so the condensation is the only possible mean to 
remove the infrared divergencies. Keeping in mind that the free energy obtained above is the finite- 
temperature part of the effective action, the appearance of the infrared divergencies at m 2 — > for 
a slowly varying gravitational field is quite expectable (see, e.g., 31-f34|). For the rapidly 



changing metric field such divergencies may not appear ||, but this fact has nothing in common 
with the infrared divergencies arising in the infrared limit. 

The coefficient of the logarithmic term is closely related with the conformal (trace) anomaly 
jl], Q. Substituting the expressions for a!f from Appendix Jb[ we come to Q [2], || [T(| 

4±^,2ft(2) ± 3 s (^ , o ~(0) , ~(2) , 3~(4) , 15~(6) , 105 ~ (8) = 



= m 4 + \m 2 R + ±(R^R^ upa - R^R^ + 6V 2 ii) + ±R 2 



A great number of cancelations happens in calculating this expression. This fact can be considered 
as the indirect verification of the correctness of the coefficients entering the expansion above. Re- 
markably, this expression does not depend on the Killing vector field. The nonminimal coupling 
adds to the anomaly (|52] ) the term, 

- 2m 2 kR + (k 2 - ^)R 2 - ^V 2 R, (53) 
3 6 

as one can easily deduce from the substitution rule (|B9[) . In the conformal case m 2 = 0, k = 1/6, 
the coefficient of the logarithm coincides with the standard expression for the conformal anomaly. 
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Some comments on the conformal transformations are in order to caution the interested reader. 
In the paper the high-temperature expansion of the free energy of a scalar field on a static 
gravitational background was obtained by the use of the conformal transform from the ultrastatic 
spacetime. Despite the fact that the final result (Eq. (18), for the terms at /3 -4 , (3~ 2 , , and 
In f3 is correct for the static gravitational field, the method for its derivation seems to be invalid 
at the point where the scaling property of the generalized zeta function (the equation before Eq. 
(10) of j|]) was employed. This scaling property holds only for the conformally invariant operator 
(m 2 = and kjj = (D — 2)/4(D — 1)) and does not take place in a general case (see, e.g., [Q). Of 
course, one can always make a transformation 

9/J.v ^ 9 fiv = 6 ^9fJ,vi 

\g\ 1 ^V 2 \g\- 1 ^- K R + m 2 = (54) 



e 



b| 1/4 V 2 M~ 1/4 — kR + {D — 1)(« - K D )(2V 2 p + (D — 2)V M/ oVV) + m 2 e 2p e~ p , 
with p independent of time. Nevertheless, the variation of this expression with respect to p gives 
not only the term arising in the scaling property of the zeta function 0], but also the contribution 
from the variation of the expression standing in the square brackets (|54|). This contribution is not 
trivial in the case k/ko and m 2 ^ 0. 

On the other hand [H], one can indeed make the substitution (|54j) from the beginning (|j) and 



apply the whole machinery of the heat kernel expansion (17) to the operator standing in the 
square brackets. The spectrum associated with this operator is obviously the same as the spectrum 
associated with the initial operator, and so the high-temperature expansions must coincide. If the 
function p is of a general type, such that, in particular, the coefficients at ui 2 and m 2 in H(uj) are 
not constants, then this function disappears from the high-temperature expansion upon restoration 
of the dependence g^ v on p. However, for a certain particular choice of p yielding a vanishing of 
some terms in the coefficients a&(o;, sc) of the heat kernel expansion, for example, for p = In 
leading to the constant coefficient at to 2 in H(uj), to the optical metric, and to the needlessness of 
the resummation of the expansion for a static metric 0, M, the dependence of the free energy on this 
p does not disappear, in general. In making the conformal transformations, one should take into 
account that the coefficients ai , j > 2, become nonvanishing and they also contribute to the free 
energy or, in other words, the additional terms arise in the scaling relation for the generalized zeta 
function |35|| . This, of course, is just a manifestation of the fact that the Klein-Gordon equation is 
not conformally invariant for the arbitrary k and m 2 . 

Now we return to the discussion of the properties of the expansions ( fi4| ) and ( ]45| ) . The contri- 
bution to the finite part of the free energy at /3 — > 0, which determines, in particular, the Casimir 
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force at zero temperature, does depend on the Killing vector. The whole expression is rather huge, 
and we do not write it here (see Eq. (B21)), but the terms at m 2 take a compact form: 



m 4 - rhHhp + jjfi< 4 >) + . . . = m 4 + ^-(R + ^R^f + ^h 2 ) + ... (55) 

These terms give the leading contribution to the finite part when m~ 1 is much smaller than the 
characteristic scale of variations of the gravitational field. The finite part does not coincide with 
the answers presented in [Q] and |^ (compare the coefficient at h 2 ). Recalling that the masses of 
all the massive particles are generated by the Higgs mechanism, we see that the term at m? gives 
a contribution to the effective potential of the Higgs field. This results in the additional correction 
to the gravitational mass-shift effect discussed in |T^| . 

If one regards the expression ( |45"|) as the regularized version of the one-loop contribution to the 
effective action at zero temperature with f3 playing the role of the regularization parameter (see 
Eq. (§)), then an interesting relation between the conformal and energy-time anomalies can be 
observed. Recall that the corresponding global symmetries lead to the (formal) Noether theorem 

s "" {x) s£ax) - 2 ^ D ". w^k) = v ^ v " r "' (56) 

where T[g^ v , 3>] is the renormalized effective action, and and are the corresponding 
currents. Usually, these relations are broken by the anomaly terms resulting from the quantum 
corrections. In our particular case (the massless conformal scalar field), we can renormalize the 
regularized contribution ( |45| ) to the effective action in several ways. 

The first one: we add the following counterterms to the initial Lagrangian 

/rW) 2 , r 2 / 2 , [- 4 + - 2 (4 2) + + ■■■] ln(/rV )• (57) 

As for the last counterterm, such a choice is unnatural and would be prohibited if the quantum the- 
ory of the vector field will be renormalizable (cf . the Coleman- Weinberg potential f36|| containing 
the term (p 4 In (ft 2 ). Nevertheless, if we add such counterterms, the global conformal invariance of the 
renormalized effective action is broken as long as this logarithmic counterterm is not conformally 
invariant. The conformal anomaly is given by the factor at the logarithm. On the other hand, the 
second relation in ( |56| ) holds on a quantum level provided the power-like divergencies are completely 
canceled out. The finite part and the coefficient of the logarithm are invariant under the stretching 
of the Killing vector. 

The second way: we add to the initial Lagrangian the counterterms 

/3" W) 2 , r 2 / 2 , [rh 4 + m 2 (4 2) + + ■■■] ln(/r 2 ), (58) 
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the coefficients at the power-like divergencies being renormalized to finite constants (not necessarily 
to zero). Then the second relation in ( |56|) is broken, while the first one is left untouched ||l5| . All the 
counterterms are invariant even under the local conformal transformations up to total derivatives. 



So, we face with a choice what the symmetry we want to preserve on a quantum level. Whatever 
the case, the dependence of the effective action on the Killing vector cannot be completely removed. 
The imaginary part of the effective action, which is responsible for the Hawking particle production 
], depends on the Killing vector. 



The field plays a distinguished role in defining the state of quantum fields. The structure 
of the counterterms suggests the simplest form of the classical action for this field. Its Lagrangian 



density reads |2||TJ, g8|, |39| 



i/ 2 -I(a|<^+ ,< + ,,■/?),/-_,/. - > (». (-,<)) 



where <f) is the Higgs field (SU(2) doublet), a, 7, and n are the dimensionless constants, and (3 



has the mass dimension 2. This action resembles the simplest model of a ferromagnetic [40| near 
the Curie point in the Landau theory of phase transitions. The constants a and /3 are such that 
the potential possesses the minimum at g 2 > 0. It is natural to suppose that the characteristic 
magnitude of the vacuum expectation value of the field g^ is of the order of the Planck scale. From 
the renormalization theory viewpoint, it is also appealing to suppose that ag 2 = /J?, where [i 2 is the 



parameter of the Higgs potential (see |12| for the notation). In this case, the constant /3 can be put 
to zero and all the coupling constants of the model become dimensionless, apart from the Newton 
and the cosmological constants. In particular, we have the estimations for the constants determining 
the potential of the field g^: a ~ r] 2 /m 2 pi and 7 ~ a 2 , where rj is the vacuum expectation value of 
the Higgs field. It should be noted that the dynamics of the field g^ are determined not from the 



least action principle, but from the self-consistency requirement [12| imposed on the effective action. 
The field g^ is at the minimum of the effective potential only for the Minkowski spacetime (see |l2|| ) . 
The exploration of the viability of the above model and its phenomenological consequences will be 
given elsewhere. Also note that the fields with higher spins can induce additional terms in (|59|). 
The spinor fields, for example, add the contribution proportional to h 2 to the terms standing in 



round brackets in (59) 
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VIII. CONCLUSION 



In this paper, we derived the high-temperature expansion of the free energy in a complete form 



). The main results of the paper are collected and discussed in Sec. |VII[ In particular, the 
explicit expressions for the divergent and finite parts of the high-temperature expansions are given 
in Eqs. (|44|), ( |45| ) in the case of the four-dimensional spacetime. We generalized the results of the 
papers [||, concerning the high-temperature expansion of the free energy on stationary spacetimes 
and found the additional contributions to the high-temperature expansion that were disregarded in 
those papers. In that part of the expansion which was presented in |9], [TIJ our expression for the free 
energy coincides with given in ||, up to integrals of the total derivatives. The method used in 
the present paper is straightforwardly generalized to the case of the higher spin fields, to the fields 
with the dispersion law differing from the relativistic one, and to the grand canonical ensembles. 



Appendix A: Useful relations with the Killing vector 

In this appendix, we collect some formulas regarding the differential calculus on the Riemannian 
manifold with the Killing vector. Let us introduce the notation: 



9fi = |f ' K = d v ln V £ 2 > fltv = (Al) 
where £^ is the time- like Killing vector of the metric g^ v . Then the following useful relations hold 
/ m a5 A = 0, g x h x = 0, = ^f^-h^g^, 

V^M" = —f^h" + g^h 2 , g^g'V^K = g 2 h 2 , (A2) 
5 A V A = -g [lM f u]x h\ V^fvx9 X = \fl„ + 9»Uxh X , V a /a^ = -\f, 
and for the curvature 



g x R\ua(i = -j^vfa(i - 2^[<x/m]^ + h^g^K + fcrfiK - g^^K, 
g x Rx» = Uxh x - \v x fx» - g^ x h\ R^g v = \f - g 2 Vxh x , 



(A3) 



where f : = = - Tr f. 

In the course of the derivation of the high-temperature expansion, it was useful to introduce the 
negative definite metric G^ u (see Eq. (|l6|)). The corresponding Christoffel symbols are 



T% = 1%, + ^(A- 1 - e)9^fl - f(h ( ,g v) g p - g^h") =: Y% + 7 £„, (A4) 
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where are the Christoffel symbols for the metric g pv . Hereinafter all the quantities associated 
to the metric G pv are designated by tildas. From (A4) we see that the connection T^ u tends to the 



projected connection (see, e.g., H, |10|, 14, |dj) when A goes to infinity. The curvature tensor for the 
connection Tf^ is defined in the standard way 

= R P ua, + + 7^7^' ( A5 ) 

After tedious calculations, we arrive at 



^h[ fl f a ] l/ g p --i 



£ 2 (9[tiK]hug f '-g[^h a] h p 'gu+^h [pj f a]u g p '-^g[,j,fa]uh p '+g[ p y a] h u g p -g [p X7 a] h p g u +f ap (h v g p -h p g u )+ 



+ f v p h [a9lA - -h [tx f a p g u + -fl^jXJ + ^(A- 1 - n 2 gug w f,ff x p - ^(A - neg^hMfgv. 

(A6) 

The Ricci curvature becomes 

iV = R ^ + li^ 1 - f)ltiv - 9(^ p f P u) - g^fu) P h p ] - eWgvih 2 - V p h p ) + g (lM f v)p h p }+ 

+ ^(A" 1 - efg^guf + + V M /i„. (A7) 

Consequently, the scalar curvature of the metric G pu is written as 

R = G pv R pv = R + 2V p h p - ^(A" 1 - e)f- (A8) 

Appendix B: Relevant parts of the heat kernel expansion coefficients 

Here we provide the parts of the heat kernel expansion coefficients relevant to the high- 
temperature expansion of the one-loop free energy in a three dimensional space. It is these 
coefficients that determine the divergent and finite parts of the high-temperature expansion when 



f3 — V 0. We just employ the known expressions |2C, 22-24] for the heat kernel coefficients adapted 



to the operator Hq (|l5|). Notice that we did not check the general formulas for the coefficients 



presented in p0| . 22 24 1 save the terms at small k. 

As follows from Eqs. (|33| ) and (^J), the relevant parts of the heat kernel coefficients au satisfy 
the inequality j > k — 4 (for d = 3), where j is the power of u entering the expression. In our case, 



using the notation of P0| , we should put 

E -»• V - E, n M „ -> -iuf^, 



^ (Bl) 
V = -K A - -V p h„, E = — -m 2 . 



l,o 1„„, <S 
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So, E has the second power in u), the curvature of the gauge fields f2 is of the first power, and the 
other structures are of the zeroth order in oj. Besides, as we resummed the heat kernel expansion 
(see Eq. (|18|)) and collected the terms (V — E) and R/6 without derivatives to the exponent, such 
terms in must be obliterated. We should to warn the reader that the Riemann tensor used in 



1 20] differs from our one by the sign, while the Ricci tensors coincide. 
Denoting the relevant part of the coefficient by tilde, we have 



d = 1, 5i = 0, (B2) 
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^R flp R tll/ + ^ -R a i3nvR al3 ^ u + |^ 2 y — \^ 2 E + gfi 2 



(B3) 



where 



V 2 R = (V 2 - hPV p )[R + 2V^ - i(A- x - £ 2 )/ 2 ], 

R^R^ = R^R^ + 2R^[h a .h v + V^K + ^(A- 1 - £ 2 )/ 2 J + 3 (A" 1 - £ 2 ) 2 [Tr f + \(f 2 ) 2 }- 
+ KA- 1 - £ 2 )[(V\/» 2 + hf 2 h + 2VW/ 2 J + (A- 1 + £ 2 )V/ PM ,F%+ 
+^/ 2 (V^ M - h 2 ) + h^V^h 2 + 2/i 2 V^ + VphuVh" - (V^) 2 , 
R^R^pa = R^R^pa - KA- 1 - eW^fpvU + (A" 1 - £ 2 )(V P / MCT V^ + 2hf 2 h)- 
-2£ 2 (VW/jL - hf 2 h + |^V M / 2 + 3h 2 f 2 + 2g 2 V„h l/ Wh» + 2g 2 h»V ll h 2 )+ 
+ |(A- 1 -a 2 [3(/ 2 ) 2 + 5Tr/ 4 ], 
V 2 V = (V 2 - hPV p )(\h 2 - ±v%)> 
V 2 E = oj 2 g 2 (6h 2 - 2V»h L 



n 2 = n^n 



-u 2 f 2 . 



(B4) 



In the expression (|B|) for d 2 and further for d k (Eqs. ©, ( ^ ), ( |B17|) , and ( |B20|) ), the indices 
are raised and lowered by the use of the metric G pv . Collecting all the above terms, we obtain 



a 2 



l 

180 



R^ pa Ruv P o + *4-R pv ^f P vf^ - R» u (Rpv + 2KK + 2VuK - Z 2 fl) + 6(V 2 - WAR 



+fe ( 5 (/ 2 ) 2 + 6Tr / 4 ) + ^(v^Uu) 2 - evpf» P ^r p - ey p fppP u K + ^nph - e^^f j 



v f2 
pu 



+ 5 4-f(Wh, 



2 J V v <"H '"J v <"H 

+WV^ v h v + M!(^V M / 2 + V 2 / 2 ) 



2/i 2 V^/i M + (V^) 2 + 10V p h u VPh u - fh^V^h 2 + lbh^V 2 K - 3V 2 V^ 



+ uj 2 g 2 
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(B5) 



The relevant part of the next heat kernel expansion coefficient is (see also [B3L |24|) 



a 3 - 3! 



^v^v^p + ^(v^) 2 + io^v 2 o^ + ±R a ^ u n a Pw + ^R^n%- 



^(V 2 ) 2 £ - V/V"£ - ^R^V^E - IVPRVpE + iV„BV"£ 



(B6) 
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where 



V^V^"' = -^[V^pVTP + 2hfh + i(A- x - £ 2 ) Tr/ 4 ], 
(V-^) 2 = - W 2 [(W^) 2 + i^ 1 " e 2 )(/ 2 ) 2 + Wf^fWh p - hfh], 
WV 2 n^ = -0J 2 [f^V 2 f pu - i(A- x - £ 2 ) Tr f - 2hfh - \hf>V p f 2 }, 
Rvp^SPPor = -u> 2 {R a ^ v f a Pr - - em 2 ? + Tr/ 4 )], 

(V 2 ) 2 £ = -2u: 2 g 2 [V 2 V^h p - &7 li h v W l h v - 6/i^V 2 /i M + 15WX7 p h 2 - 5h^V p V u h u - 2(V^ - 3/i 2 ) 2 ] 
i^ft 2 , = -u^BTft, + KA- 1 - £ 2 ) Tr / 4 + /i/ 2 /* + f 2 u Wh% 
V„VWE = u 2 g 2 [WV p V u h» - \WV„h 2 \, 

R^WV^E = uj 2 g 2 [R pu (4h»h u - 2V»h v ) + (3A~ X - f)hfh + 6/i 4 + h»V p h 2 - 2V p h u V»h u - 

-2h 2 v^ - (a- 1 - exfawh" +v»f ia ,r<>h p ) + \h 2 ef\ 

WRV p E = -2u 2 g 2 h v V v [R + 2Wh„ - ^(A" 1 - £ 2 )/ 2 ], 
V p EV^E = Auj 4 g 4 h 2 . 

(B7) 



a 3 = J - 



And so 

_ [^-R^f^f^ + fl^/jL + - 2VM + 2f (Tr f + |(/ 2 ) 2 ) 

+2£ 2 V^Wr" + f(V7„) 2 + 3£ 2 /^V 2 / M „ - ^/i/ 2 /* + 2£ 2 Vf U(1 f»Ph p + 3£ 2 /^V,J 2 
+2£ 2 V< 1 h u f 2 u - 3V 2 Vh p + WV»h u X7 p h u + 18/^V 2 /i M - f /i^V^/i 2 + 18h^V fl V u h u + 6(W/i M ) 2 
-38/i 2 V^/i M + 60/i 4 - ^f/i 2 / 2 - 6h^V p R 



+ I^V/i 2 . 

(B8) 

If we had included the nonminimal coupling kR of the scalar field with gravity to the Klein-Gordon 

equation (Q), this just would have resulted in a change of the coefficients of the scalar curvature 

entering the expressions for d 2 and 03 and a redefinition of the effective mass rh: 

2 2 

rh 2 ^m 2 - kR, a 2 -> a 2 - -(V 2 - h^Vn)R, a 3 -> a 3 - - 9 K h>*V u R. (B9) 

6 3 

The higher coefficients are encrypted in [ ^| . Notice that formula (6.6) for Z( 2 ) in [22| contains a 
spurious term i? 2 /12. Making use of the notation from that paper, we write out the relevant parts 
of the coefficients: 



04 — 3f 



t4iZZi) + 2 { zf lZ3) } + |zt 2 Z 2) + IZ 2 2) , (BIO) 



where 



^(l^ 1 ) ~ ~~ \R^ V x p x u + Ix^fi^O"^, 

{Z^} « + §x"x^ + A^(tt 2 ) ; „ - f x^n^ A , 

Z{ 2 Z 2) « ±[(A^) 2 + 2X^J^ - 2A^ 2 „ + X» p n 2 + i(0 2 ) 2 + iTr0 4 ], 

z 2 2) -(n) 2 +n^ 2 +i(^ 2 ) 2 - 
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Here X = V — E as above, the indices of X and semicolon denote the covariant differentiation with 
respect to the Levi-Civita connection Tfu,, and the approximate equality means that the relevant 
part of the expression is retained. Therefore we have (see also [[24||) 



a 4 — IT 



where 



1 V 2 EQ 2 + ±(n 2 ) 2 + AWVEVuVvE + &WV v ESl%, + 4 Tr ft 



>iv i 15 



(B12) 



R^V^EV^E = 2u A g A [2R flu h> 1 h u + (A" 1 - £ 2 )/i/ 2 /i + 2/i 4 + /i^V M /i 2 ], 

v^n^n^VpS = 2 W y (V7^r^ p - /i/ 2 /!), 

WEVf&E = 4^V(/i' t V M V ! ,/^ - 2h 2 V tI h ll + 6h A - 3/i^V^ 2 ), 

V^EV^n 2 = 2wV^V /i / 2 , 

(V 2 E) 2 = 4^V(3/i 2 - V^) 2 , 

V 2 £ft 2 = 2wV/ 2 (V M ^ - 3/i 2 ), 
(n a )2 = w 4 (/2)2) 

V^V u EV tl V v E = 4w 4 # 4 [V^W/i^ - 2/^V M /i 2 + 3/i 4 - ^(A _1 - £ 2 )hf 2 h], 
WV^EV 2 ^ = 2co A g 2 (Wh»fi iy ~ 2h f 2fl )> 
Trfi 4 =w 4 Tr/ 4 . 



(B13) 



Summing up, we obtain 
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iV W + U^f'Y + 4Tr / 4 ) - 7 -^hf 2 h + ^Wfc + " T V^Uuf up h p - ^Vj 2 

lf/ 2 (V^ - 3h 2 ) + fh 4 - fWV^h 2 + 2WV^ u h» - 9/i 2 V^ t + §(V%) 2 + f V^VW] . 

(B14) 



Further, 



°5 - kt 



S Z iZ 2z D + WMi) + UAi 1 ^ + H z fi Z D Z W> 



(B15) 



where 



Z fi Z2Z i) 



x^ 2 „x> 



Z^ZaZjj « fX^X^X^ + \X^X V V + ±X^X M ft 2 - IX^X", 



{z| 1 ZiZ 2) }^2Zj 1 Z 2 Z 1) , 

{Z^ZxjZp)} « 2[x^x M x j ; + ix^x M o 2 

Hence, substituting ( gig ) to ( pl5| ), we come to 



(B16) 



a 5 - 5, 



2V^EV U EV.V U E - \V^EVaEV 2 E + ^V^EV.En 2 - \V^EV v E^l 



(B17) 
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where 



V^EV U EV ^V V E = 4w 6 c/ 6 (4/i 4 - /i^V^/i 2 ), 
Vt*EV^EV 2 E = 8uj 6 g 6 h 2 (3h 2 - V/» M ), 
V^EV^EQ 2 = -4u 6 g 4 h 2 f 2 , 
V f "EV u Enf lu = -4uj 6 g 4 hf 2 h. 



Collecting all the terms together, we get 



(B18) 



a 5 



4wV 



18/i 4 - 2/i^/i 2 - ftfVhn + 5 -fh 2 f 2 - ^hfh 



(B19) 



The last coefficient that we need is 



«6 — el 



jX^X^X^X^ + j^X^XiXiX^ 



±l(WEV ll E) 2 = §u 8 g 8 h±. 



(B20) 



In conclusion we give the explicit expression for the finite part of the high-temperature expansion 
of the one- loop contribution to the free energy (see Eqs. (p4|), (45)) retaining only the terms at the 
nonegative powers of m 2 : 



^4 _ m 2 (|4 2) + |4 4) ) - |(4 2) + 2~af ! + f~a ( °> + 225 



;( 4 ) , 23- (6) 



(8)v 



m 



+ 



15 
.4£ 



4 + m 2 {\R+%f 2 -^h, + %h 2 ) 

e^^fapf^ + w v (ef% - - 2h,h u ) + ^(Tr/ 4 + §(/ 2 ) 2 ) + ^(v^/^) 2 

+^F(v M /, p v"r p - - /£,vw) + 2e(P u v 2 Uu + fc"V M / 2 ) - §/x 2 / 2 + ±fV 2 v% 

-2V a Vfc„ + fWWVuK + 8/i^V 2 /i M + ±(V^ t ) 2 + f /i 2 V^V - f/^V^/i 2 
" i^V^ + iflf ^ 2 - + f f 2 ) + ±R 2 . 
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The nonminimal coupling adds to this part the following terms: 

_ K R(^h 2 - -V»K + —f + 2m 2 ) + {k 2 - -)R 2 + —h^V.R. 
48 9 " 36 J ' v 3 ; 9 M 



(B21) 



(B22) 



In the conformal case, m 2 = 0, K = 1/6, this contribution cancels completely the terms proportional 
to the scalar curvature -R (without derivatives) entering the finite part. 
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